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The semiclassical, dynamical description of diffraction dissociation of weakly bound nuclei is 
applied to analyze the decay-energy spectra of 15 C that have been measured at 68 MeV/nucleon on 
a Pb target. The optical potentials that are used to describe the nuclear interaction of 15 C with 
the target nucleus are realistic because the fits to the two measured spectra, one with a small and 
one with a very large acceptance angle, are consistent and of similar quality. The cross section for 
the radiative neutron capture on 14 C to the l/2 + ground state of 15 C is deduced from the analysis. 
When combined with an estimated contribution from the capture to the 5/2 + excited state of C, 
an excellent agreement with a recent direct capture measurement is achieved. 

PACS numbers: PACS number(s): 25.40.Lw, 25.60.Gc, 25.70.De, 27.20.+n 



I. INTRODUCTION 

The decay-energy spectrum for the breakup of 15 C on 
a Pb target was first reported in Ref. The data 

were recently analyzed by CDCC (continuum discretized 
coupled-channels) calculations 0, and it was demon- 
strated that the measured spectrum is consistent with 
measurements of the radiative neutron capture rate on 
14 C [3j. The final results of the Coulomb dissociation ex- 
periment have now been published (4j. They include the 
decay-energy spectrum for all events and a spectrum for 
events at forward angles. 

It is of interest to investigate whether the two spectra 
obtained in Ref. [i[ can be reproduced by a single, con- 
sistent theoretical description, because the forward angle 
measurement is primarily due to Coulomb dissociation, 
whereas events at large scattering angles are influenced 
by nuclear processes. A more general issue is whether 
the cross section for the radiative neutron capture on 
14 C that can be deduced from the Coulomb dissociation 
experiment is consistent with the direct capture measure- 
ment of Ref. 3 . 

The method that will be used to analyze the decay- 
energy spectra of 15 C is the semiclassical, dynamical 
description, where the relative motion of the projectile 
and target follows a classical Coulomb trajectory. The 
breakup reaction, 15 C— > 14 C+n, is calculated quantum 
mechanically by solving the time-dependent Schrodinger 
equation for the relative motion of the valence nucleon 
and the core fragment, in the time-dependent Coulomb 
and nuclear fields of the target nucleus. 

The semiclassical model used here was first introduced 
in Ref. Q as a three-dimensional generalization of a pre- 
vious two-dimensional model [6]. It has been applied in 
studies of the Coulomb dissociation of 8 B 0,0 and 17 F@ 
nuclei, and it has provided a qualitative understanding 
of some of the phenomena that can occur. For example, 
intermediate energy Coulomb dissociation experiments 
are commonly analyzed in first-order perturbation the- 
ory in terms of El transitions using the so-called far-field 
approximation and straight line trajectories [ioj . These 
approximations are rather poor for the Coulomb disasso- 



ciation of weakly bound proton+core systems, partly be- 
cause the far-field approximation breaks down and partly 
because E2 transitions and higher-order processes can- 
not be ignored SI- As a consequence, a discrepancy of 
15% was observed between the measured radiative cap- 
ture rate of protons on 7 Be and the rate inferred from 
Coulomb dissociation experiments [ll[ . The discrepancy 
was later resolved by CDCC calculations [12] that in- 
cluded the effects discussed above, in addition to the nu- 
clear induced breakup of 8 B. 

It is surprising that the discrepancy between the di- 
rect measurement of the cross section for the radiative 
neutron capture on 14 C Q and the first-order analysis 
of the 15 C Coulomb dissociation experiment Q is also 
of the order of 15%. Naively, one would expect that the 
analysis of the measured decay-energy spectra in terms of 
first-order perturbation theory Q would result in a much 
better agreement with the direct capture measurement 
because the far-field approximation for Coulomb excita- 
tion is well justified and the E2 strength is very weak 
for a neutron halo nucleus. Consequently, the dynamic 
polarization effect, which is caused by an interplay of El 
and E2 transitions [13], is much weaker for a neutron 
than for a proton halo nucleus. It is therefore of interest 
to analyze the measured decay-energy within the semi- 
classical, dynamical description, in order to see whether 
the discrepancy with the radiative capture measurement 
can be reduced. 

Part of the explanation for the 15% discrepancy is that 
the measured neutron capture rate, in addition to the 
direct capture to the 1/2+ ground state of 15 C, also in- 
cludes a small branch to the 5/2 + excited state, whereas 
the Coulomb dissociation experiment can only provide in- 
formation about the capture to the ground state. Unfor- 
tunately, the cross sections for the capture to the ground 
state and the excited state have not been measured sep- 
arately. In this work the contribution from the 5/2 + 
branch is estimated to be of the order of 4%, so the issue 
is what causes the remaining 11% discrepancy. 

Details of the structure input and some results of first- 
order perturbation theory are presented in the next sec- 
tion. The basic ingredients of the semiclassical method 
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TABLE I: The single-particle structure of 15 C is simulated by 
a Woods-Saxon (WS) potential. The radius R and diffuseness 
a are fixed, and the depth is adjusted for each partial wave, 
V s , V p and Vd, so that certain properties of 15 C are reproduced 
(see the text.) Also shown is the spin-orbit strength, V so . 



R (fm) 


a (fm) 


V 3 (MeV) 


V p (MeV) 


V d (MeV) 


Vso (MeV) 


2.946 


0.5 


55.36 


55.36 
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FIG. 1: (Color online) The calculated dipole and quadrupole 
responses of 15 C are shown as functions of the relative energy 
Erei of the neutron and the 14 C fragment in the final state. 



are summarized in Sect. III. The convergence of the cal- 
culations and the analysis of the new decay-energy spec- 
tra is presented in Sect. IV. The results are used in Sect. 
V to infer the cross section for the radiative neutron cap- 
ture on 14 C. The conclusions are presented in Sect. VI. 
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FIG. 2: (Color online) Coulomb dissociation cross sections 
for 15 C on a Pb target as functions of the beam energy. The 
results of first-order dipole excitations are shown for straight- 
line, corrected straight-line (& e // = b + Tia/2), and for exact 
Coulomb trajectories. Dipole excitations up to E re i = 4 MeV 
are included, and the distance of closest approach is larger 
than Rx+R 2 = 1.2(A 1/3 + A\ /3 ). 

No choice was made in Ref. [|| for the d-wave potential 
because it was not needed in the first-order calculations of 
dipole excitations performed there. The depth shown in 
Table I has been determined in this work so that the 0.478 
MeV binding energy of the 5/2 + excited state of 15 C is 
reproduced. Combined with the spin-orbit interaction, 
this implies that the 3/2 + resonance is located in close 
vicinity of the known resonance at 3.56 MeV. 

A. First-order Coulomb excitation 



II. STRUCTURE MODELS OF lj C 

The single-particle structure associated with the va- 
lence neutron in 15 C is simulated by a Woods-Saxon 
(WS) potential. The depth is adjusted for each partial 
wave so that certain properties of 15 C are reproduced. 
The parameters of the model are shown in Table I. The 
radius and diffuseness, i?=2.946 fm and a=0.5 fm, were 
actually determined in a previous work [l4| by simulat- 
ing the mean-field potential of a realistic Hartree-Fock 
calculation. The depth of the s-wave potential, V s , was 
adjusted to reproduce the 1.218 MeV neutron separation 
energy of the l/2 + ground state. 

The depth of the p-wave potential, V p , is assumed to be 
the same as the depth of the s-wave potential. This choice 
was made deliberately in Ref. Q , in order not to violate 
Siegert's theorem which allows one to use charge densities 
instead of current densities when calculating electromag- 
netic matrix elements (see, e. g., Ref. J15j). The depth of 
the WS potential for higher, odd-parity partial waves is 
chosen in the following to be the same as for the p-wave 
potential. 



The El and E2 strength distributions for 15 C that one 
obtains are illustrated in Fig. Q] as functions of the fi- 
nal state relative energy E re i of the neutron and the 14 C 
fragment. The E2 response has a shoulder between 3 
and 4 MeV which is a consequence of the 3/2 + reso- 
nance. The E2 response does not play any practical role 
in the Coulomb dissociation of 15 C because the E2 exci- 
tation probability is more than three orders of magnitude 
smaller than the El excitation probability. However, the 
3/2 + resonance does play a role in the nuclear induced 
breakup. 

Coulomb dissociation experiments at intermediate and 
high energies are usually analyzed using the first-order 
perturbation theory that is based on straight-line trajec- 
tories [10j. A simple correction for Coulomb trajectories 
is often made in this description. It consists of replac- 
ing the minimum impact parameter bo for a given mea- 
surement by the effective minimum impact parameter 
bo + na/2, where a is half the collision diameter. This is 
actually a very good approximation for first-order dipole 
excitations. This can be seen in Fig. [5] by comparing 
the middle dashed curve (labeled b e ff = b + na/2) to the 
solid curve, which is the first-order Coulomb dissociation 
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cross section obtained using exact Coulomb trajectories 
as described in Ref. [HI]. 

The calculations shown in Fig. [2] were based on the 
dipole response shown in Fig. [1] and included relative 
energies up to 4 MeV and all impact parameters for 
which the distance of closest approach is larger than the 
sum of the radii, J? x + R 2 = 1.2(A* /3 + A 1 / 3 ). The top 
dashed curve is the cross section for straight-line trajec- 
tories. This is a poor approximation for energies below 60 
MeV/nucleon. The corrected straight-line trajectory ap- 
proximation (middle curve) is in much better agreement 
with the calculation that is based on exact Coulomb tra- 
jectories (the solid curve). At 68 MeV/nucleon, the two 
calculations agree within 1%, whereas the pure straight- 
line trajectory calculation (the upper dashed curve) is 
about 5% higher. These results show that the corrected 
straight-line trajectory approximation, which was used in 
the experimental analysis is well-justified, provided 
higher-order and nuclear processes can be ignored. 

III. THE SEMICLASSICAL METHOD 

The relative motion of projectile and target is assumed 
to follow a classical Coulomb trajectory in the semiclas- 
sical description of breakup reactions. The trajectory in- 
cludes relativistic effects, the most important being the 
kinematics and the determination of the velocity from 
the beam energy. These aspects are discussed in detail 
in Ref. [H]. 

The relative motion of the valence neutron and the 
14 C core is calculated by solving the time-dependent 
Schrodinger equation [5(. The two-body wave function 
^(r,t), which has the l/2 + ground state as initial condi- 
tion, is expanded on the angular momentum eigen states 
\ljm), 

1 L max 

*( r >*) = ~ 12 u um(r) l*i m >- (!) 

Ijm 

The upper limit L max is the maximum orbital angular 
momentum of the expansion. The radial wave functions 
uij m (r) are calculated on a radial grid out to 100 fm, with 
a grid size of 0.1 fm. The total wave function is evolved 
in time by the propagator 

[l-l^o]" 1 [l+^o] *(*), 

(2) 

where Ho is the two-body Hamiltonian for 15 C, and 
V ex t(t) is the time-dependent interaction of 15 C with the 
target nucleus. The form of Eq. (2) was originally used 
m Ref. 0. The propagator that was used in the earlier 
work [f| was a simplified version of Eq. (2). 

Inserting the expansion (1) into Eq. (2) one obtains 
the following expression for the propagation of the radial 
wave functions, 



ui jm (r,t + St) = [l- 




(uij m (r,t) + — 12 {ljm\Vext\l'frn') ui'j/ m >(r,t)), (3) 

I'j'm' 

where hij is the radial single-particle Hamiltonian which 
depends on (Ij) through the spin-orbit interaction. The 
coupling between the radial wave functions occurs only 
through the interaction with the target, whereas the uni- 
tary, intrinsic propagation is diagonal in \ljm). The in- 
verse operator, [l — St/(2iK)hij]~ 1 , is calculated using the 
technique described in Appendix B of Ref. [IH . 

The spin-orbit interaction is included explicitly in the 
following. That makes it possible to treat certain aspects 
of 15 C in a realistic way, such as the 5/2 + bound excited 
state and the 3/2 + resonance. The interaction V ext be- 
tween the projectile and target consists of the core-target 
Coulomb interaction, the core-target nuclear interaction, 
which is based on the 17 0+Pb optical potential of Ref. 
[l9l ] , and the Perey-Perey neutron-target interaction [2(| • 
These interactions were chosen here because they were 
applied in the CDCC calculations of Ref. @. They were 
expanded on Legendre polynomials, 

V ext (r,t) = 12 Vext,x(r,t) Pa(cos(0')), (4) 

A=0 

where 0' is the angle between the position r of the neu- 
tron with respect to the 14 C core and the trajectory R(i) 
of the target with respect to the 15 C projectile. The max- 
imum value of A in the expansion \ m ax, was set equal 
to 2L max so that all of the necessary multipole compo- 
nents were included in the calculation of the matrix ele- 
ments of V ex t that appear in Eq. ([3]). 

The semiclassical method is applied to the breakup of 
15 C on a Pb target at a beam energy of 68 MeV/nucleon. 
The time-evolution of the wave function for a given im- 
pact parameter b starts when the distance between pro- 
jectile and target is 300 fm with the 1/2+ ground state 
wave function as the initial condition. The propagation is 
terminated for practical reasons when the projectile and 
target have re-separated by 100 fm. One could follow 
the evolution further but the result would be useless be- 
cause the wave function would start to be reflected from 
the outer boundary of the radial box and that produces 
oscillations in the calculated decay-energy spectrum. 

One way to avoid the problems caused by the reflection 
of the wave function from the outer boundary is to use a 
larger radial box. This method was used to test the con- 
vergence of the calculated decay-energy spectrum, and 
it appears to have converged already at a separation of 
100 fm. Another way of avoiding the reflection from the 
outer boundary (which will not be used here) is to apply 
an imaginary potential of a special form that acts near 
the outer boundary; see Ueda et al. [2l[ for details. 

IV. APPLICATIONS TO 15 C BREAKUP 

In this section the semiclassical method is applied to 
analyze the decay-energy spectra of 15 C that were mea- 
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FIG. 3: (Color online) The diffraction dissociation probabil- 
ity Pdiff(b) for 15 C on a Pb target at 68 MeV/nucleon (solid 
red curve) is shown as a function of the impact parameter b 
and is compared to the first-order El Coulomb dissociation 
probability (black dashed curve). The steeply falling (green) 
dashed curve shows the total absorption probability. 
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FIG. 4: (Color online) Decay energy spectra at the impact 
parameter b — 12, 20, and 30 fm. The solid curves are dy- 
namic calculations, whereas the dashed curves are first-order 
perturbation calculations. 



sured at 68 MeV/nucleon on a Pb target Q. The pa- 
rameters of the structure model are shown in Table I; 
they are essentially the same as those that were used in 
the previous first-order analysis of the measurement [3| • 
Before presenting the results of the new analysis of the 
data, it is useful first to show the dissociation probabil- 
ities one obtains, and also to study the convergence of 
the calculations with respect to the maximum angular 
momentum L max that is used in the expansion (1) of the 
15 C two-body wave function. 



A. Disassociation probabilities 

The calculated dissociation probabilities are shown in 
Fig. [3] as functions of the impact parameter. The solid 
curve is the dynamic calculation with L max =4. The 
dashed curve (labeled El Coulex) is the result of the 
first-order El Coulomb excitation, which is here based 
on Coulomb trajectories [l6|. The dynamic calculation 
is slightly suppressed at large impact parameters com- 
pared to the first-order calculation. In the 11-16 fm im- 
pact parameter range, the nuclear induced breakup takes 
over and gives an enhancement compared to the first- 
order calculation. Finally, at impact parameters smaller 
than 10 fm, the imaginary core-target interaction be- 
comes strong and absorbs most of the wave function. 
The steeply falling dashed curve is the total absorption 
probability, which is due to the combined effect of the 
imaginary part of the neutron-target and core-target in- 
teractions. 



B. Convergence of the dynamic calculations 

The decay-energy spectra obtained at three impact pa- 
rameters are compared in Fig. [4] to those obtained from 
first-order El transitions. The spectra obtained in the 
dynamic calculations are very broad at small impact pa- 
rameters but they approach the first-order calculation 
at the large impact parameters. The decay-energy spec- 
trum, da/dE, will therefore in the following be calculated 
numerically by integrating the dynamic spectra over im- 
pact parameters less than 30 fm, whereas the contribu- 
tion from impact parameters larger than 30 fm will be es- 
timated by the first-order Coulomb excitation spectrum. 

The dependence of the dissociation probability on the 
maximum angular momentum L max is illustrated in Fig. 
[5] in terms the continuum occupation probabilities P(L) 
(summed over j = L± 1/2) for a fixed impact parameter 
of 12 fm. It is seen that the occupation probability for 
a fixed value of L is largest when L max =L, but the value 
drops and converges rather quickly with increasing val- 
ues of L max . The total dissociation probability is 0.107, 
0.102, 0.097 and 0.096 for L max = 2, 3, 4, 5, i. e., the 
total probability is reduced by about 10% by increasing 
Lmax from 2 to 4. The occupation probability of the 5/2 + 
bound state is also shown in Fig. Oat L=2; it is relatively 
small and converges to a value of about 0.002. 

The decay-energy spectra one obtains in the dynamic 
calculations are illustrated in Fig. [5] for two values of 
the maximum angular momentum, L max = 2 and 4, re- 
spectively. The peak height decreases from 393 to 375 
mb/MeV as L max increases from 2 to 4; that is a reduc- 
tion of 5%. The peak height is not expected to change 
much by increasing the value of L max further because the 
reduction is only 1% when L max is increased from 3 to 4. 

The nominal maximum center-of-mass acceptance an- 
gle in Fig. [S]was set to 6° because that is the largest ac- 
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FIG. 5: (Color online) Occupation probabilities of the con- 
tinuum as function of the orbital angular momentum L of the 
neutron+ 14 C system, at the fixed impact parameter 6=12 fm. 
The probabilities are shown for different values of the maxi- 
mum orbital angular momentum, L max . The occupation prob- 
ability of the excited 5/2 + bound state (BS) is indicated for 
L—2 inside the lower left box. 
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FIG. 6: (Color online) Decay energy spectra obtained in dy- 
namic calculations with L maa; = 2 and 4, respectively. The 
nominal acceptance angle is 6° but the calculated spectra have 
already converged at a smaller angle (see text.) 



ceptance angle that was used in the experiment [J]. The 
dynamic calculations have actually already converged at 
a smaller scattering angle because of the strong absorp- 
tion at small impact parameters. For example, at the 
impact parameter b — 9 fm, where the dynamic disso- 
ciation probability is already quite small, the Coulomb 
scattering angle is 6* cm = 4.9 °. 



C. Analysis of the measured spectra 

The decay-energy spectrum of 15 C at 68 MeV/nucleon 
on a Pb target was measured for two center-of-mass ac- 
ceptance angles, namely, at 2.1° and 6° [4J. The data 
are compared to the first-order and the dynamic calcu- 
lations in Fig. and [7)3, respectively. The actual cal- 
culations are shown by dashed curves. The solid curves 
were obtained by a folding and scaling procedure: the 
calculated spectra were first folded with the experimental 
energy resolution, which is a Gaussian with a (la) width 
otAE rel = Q.23VE^i 0. The folded spectra were next 
scaled by the factor S c which optimizes the x 2 fit to the 
data. The values of the scaling factors and the associated 
best values of the \ 2 P er point are listed in Table II for 
different values of L max . 

The decay-energy spectra for 9 < 2.1° were calculated 
with a sharp cutoff at the impact parameter where the 
Coulomb scattering angle is 6 = 2.1°. An acceptance an- 
gle of 6° does not impose any strict cutoff on the dynamic 
calculation because it converges already at a smaller an- 
gle. The first-order Coulomb excitation calculation, on 
the other hand, does not have such a natural cutoff, ex- 
cept from Coulomb scattering but that leads to a very 
large cross section. The first-order decay-energy spec- 
trum (El Coulex) shown in Fig. [JJ^ for 9 < 6° was 
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FIG. 7: (Color online) Decay energy spectra obtained in first- 
order perturbation theory (A) and in the dynamic calcula- 
tions with Lmax=4: (B) are compared to the data of Ref. [j. 
The dashed curves show the calculated spectra. The solid 
curves include the experimental energy resolution and have 
been scaled to give an optimum fit to the data. 
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TABLE II: Analysis of the measured decay-energy spectra of 
15 C on a Pb target [J. The scaling factor S c and the asso- 
ciated x 2 /N for the best fit to the data up to E re i < 4 MeV 
are shown as functions of the maximum angular momentum 
Lmai, and for the two acceptance angles of the experiment. 
Also shown are the cross sections before and (in parenthesis) 
after the folding and scaling procedure. 





9 < 2.1° 
a (mb) 


9 < 2.1° 

Sc - X 2 /N 


9 < 6.0° 
a (mb) 


9 < 6.0° 
S c - X 2 /N 


El Coulex 
2 
3 
4 


326 [303] 
316 [303] 
316 [303] 
316 [304] 


0.941 - 0.87 
0.972 - 0.76 
0.972 - 0.77 
0.973 - 0.77 


767 [638] 
750 [696] 
732 [699] 
716 [696] 


0.841 - 2.24 
0.939 - 1.15 
0.969 - 1.20 
0.984 - 1.19 



therefore determined by integrating over all impact pa- 
rameters for which the minimum distance of closest ap- 
proach is larger than 1.2( J 4 1 / 3 + A^ z ). The associated 
minimum impact parameter for a Coulomb trajectory is 
9.7 fm, which is a fairly reasonable choice because it falls 
in the region where the absorption in the dynamic calcu- 
lation sets in (see Fig. [3]) 

The dynamic calculations converge quickly for the 
smaller acceptance angle. This can be seen in Table II 
where the scaling factor S c that gives the best fit to the 
data is independent of L max . At the 6° acceptance an- 
gle, the scaling factor increases by almost 5% to the value 
S c w 0.98 for L max =4; it is not expected to increase much 
further for larger values of L max . The scaling factors ob- 
tained in the two analyses are therefore approximately 
identical and the x 2 /N is also very reasonable for both 
acceptance angles. This implies that the adopted nuclear 
interactions with the target must be realistic because the 
calculations at large scattering angles are strongly influ- 
enced by the nuclear interactions, whereas the dissoci- 
ation at the smaller acceptance angle is dominated by 
Coulomb dissociation. 

The results of the first-order analysis are shown in the 
first line of Table II. The analysis of the large acceptance 
angle measurement is not so interesting because the fit 
is poor and the necessary scaling factor is small and un- 
certain. The uncertainty stems from the crude estimate 
of the minimum impact parameter. Although one could 
possibly choose a better value for the minimum impact 
parameter, the fit to the data would still be poor because 
the first-order decay-energy spectra are narrow compared 
to the results of the dynamic calculations at small impact 
parameters. This can be seen in Fig. 0J 

The result of first-order perturbation theory at the 
smaller acceptance angle is much more interesting. Here 
the x 2 '/N is good and the necessary scaling factor does 
not differ dramatically from the dynamic calculations. 
However, we shall see in the next section that the 4% 
larger scaling factor of the dynamic calculation, com- 
bined with other corrections, is essential for reaching a 
good agreement with the neutron capture data. 

It should be emphasized that the scaling factors S c 



TABLE III: MAC cross sections (kT = 23.3 keV) for the 
radiative neutron capture on 14 C to the l/2 + ground state 
of 15 C, to the 5/2 + excited state, and the sum. The first 
line is the prediction of the structure model (with S c = 1.) 
The second and third lines are the results of the first-order 
El and the dynamic calculation analysis of the decay-energy 
spectrum with a 2.1° acceptance angle. The last line is the 
measured cross section [3|. 



Method 


s c 


(tmac(1/2 + ) 


(tmac(5/2+) 


ctmac (total) 






( M b) 


(Mb) 


( M b) 


Model 


1 


6.85 


0.26 


7.11 


El Coulex 


0.941 


6.45(50) 


0.24 


6.69(50) 


Dynamic 


0.973 


6.67(50) 


0.25 


6.92(50) 


Exp. [3] 








7.1(5) 



listed in Table II should not be confused with spectro- 
scopic factors. The scaling factors were used here as a 
convenient way to analyze the data and show how the 
calculations converge with increasing values of L max . 



V. COMPARISON TO RADIATIVE CAPTURE 

As mentioned in the introduction, the discrepancy be- 
tween the direct and indirect measurements of the ra- 
diative neutron capture rate on 14 C is about 15%. This 
statement is based on the values of the Maxwellian Av- 
erage Capture (MAC) cross sections, ctmac, that were 
obtained at the temperature kT = 23.3 keV. The defini- 
tion of the MAC cross section is quoted in Eq. (Al) of 
the appendix. The direct measurement 0] gave the value 
<tmac = 7.1 ± 0.5 jjb, whereas the indirect, first-order 
Coulomb dissociation method gave the cross section 6.1 
± 0.5 /A> [H]. One reason for the discrepancy is that the 
direct measurement includes a small contribution from 
the capture to the 5/2 + excited of 15 C, whereas the in- 
direct measurement can only provide information about 
the capture to the l/2 + ground state of 15 C. 

In order to investigate in more detail what causes the 
discrepancy between the Coulomb dissociation and ra- 
diative capture measurements, it is useful to have an es- 
timate of the cross section for the radiative capture to the 
5/2 + excited state of 15 C. The prediction of the structure 
model that was used in the previous sections is shown in 
the first line of Table III. The total MAC cross section 
predicted by the model is (accidentally) in perfect agree- 
ment with the measured 7.1(5) fib cross section ||. The 
contribution from the capture to the 5/2 + excited state 
is almost 4%, so the discrepancy between the first-order 
Coulomb dissociation analysis of Ref. [1] and the radia- 
tive capture measurement is reduced to 11%. 

The second and third lines of Table III show the cross 
sections extracted from the analysis of the decay-energy 
spectrum that was measured with the 2.1° acceptance 
angle (see Table II.) The cross sections were obtained by 
multiplying the predictions of the structure model (the 
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10 100 1000 

E rel (keV) 

FIG. 8: (Color online) Calculated cross sections for the radia- 
tive neutron capture on 14 C to the l/2 + ground state (G.S.), 
the 5/2 + bound state (B.S.) of 15 C, and the sum (Total) are 
shown as functions of the 14 C+neutron relative energy, E re i. 
The measured cross sections are from Ref. Q. 

first line of Table III) with the scaling factors that gave 
the best fit to the measured spectrum. The cross sec- 
tions have been assigned an uncertainty of 0.5 /ib, which 
roughly reflects the 7% uncertainty of the Coulomb disso- 
ciation experiment The total cross section obtained 
from the first-order analysis (second line) is smaller than 
but not inconsistent with the capture measurement. The 
cross section extracted from the dynamic calculation 
(third line) is in excellent agreement with the measure- 
ment. Thus it is the combined effect of dynamic processes 
and the estimated contribution from the capture to the 
5/2+ state that makes it possible to achieve consistency 
between the Coulomb dissociation and the direct capture 
measurements. 

The MAC cross section for the radiative capture to the 
1/2+ ground state of 15 C that was obtained in Ref. [1| 
is 6.1(5) /ib. That is about 5% smaller than the 6.45(50) 
/Ltb cross section shown in the second line of Table III. 
The two results should in principle be identical because 
they are based on the same data set (the decay-energy 
spectrum at 2.1°) and they were both obtained using 
first-order perturbation theory. In this connection, it is 
very interesting that the two analyses appear to give the 
same dipole strength distribution. For example, the peak 
height of the extracted distribution shown by the dashed 
curve in Fig. 2 of Ref. [| is 0.32 e 2 fm 2 /MeV. The same 
peak height is obtained from the dipole response shown 
Fig. 1 when multiplied with the (El Coulex) scaling fac- 
tor of Table II: 0.34 x 0.941 = 0.32 e 2 fm 2 /MeV. The 5% 
discrepancy mentioned above must therefore have devel- 
oped in the calculation of the MAC cross section from the 
extracted dipole strength distribution (see appendix.) 

The radiative capture cross sections predicted by the 
structure model defined by Table I are shown in Fig. [8] 
as functions of the neutron center-of-mass energy. The 
figure shows the contributions from the capture to the 
1/2+ ground state and to the 5/2 + excited state of 15 C. 



The solid curve is the sum (total) which can be compared 
to the cross sections that were measured in Ref. @ • The 
model is in very good agreement with the data. The best 
prediction obtained from the analysis of the Coulomb dis- 
sociation experiment is a factor of 0.973 times the model 
prediction which is also in good agreement with the data. 

The experimental cross section shown at the lowest 
energy (23.3 keV) in Fig. g]is 4.86(34) /ib. This value was 
determined from the published 7.1(5) /ib MAC cross and 
the simple relation, Eq. (A3), derived in the appendix. A 
different value, namely 5.2(3) /ib, was quoted in Ref. [|[ 
but that value has been discarded here because it does 
not agree with the prediction of Eq. (A3), which was 
derived from the structure model used here and ei ghth 
other structure models that were considered in Ref. 12211 . 



VI. CONCLUSIONS 

The semiclassical, dynamical description of the disso- 
ciation of weakly bound nuclei was applied to analyze 
the decay-energy spectra of 15 C that have been measured 
in reactions with a Pb target. The two-body structure 
model of 15 C, which was partly developed previously, was 
supplemented with a model for d-waves, so that the en- 
ergies of the 5/2 + bound state and the 3/2 + resonance 
could be simulated. Standard neutron-target and core- 
target optical potentials were employed to describe the 
nuclear interaction of 15 C with the target nucleus. 

The two spectra that were measured were obtained 
with a small and a large acceptance angle, respectively. 
The calculated spectrum for the small acceptance an- 
gle is dominated by Coulomb dissociation, whereas the 
spectrum calculated with the large acceptance angle is 
strongly influenced by the nuclear interaction. The ef- 
fect of the nuclear interaction is to produce broad decay- 
energy spectra at large scattering angles, and it also de- 
termines the magnitude of the spectrum for the large ac- 
ceptance angle because of a strong absorption it produces 
at small impact parameters. 

The analysis of the measured decay-energy spectra 
shows that the optical potentials that were used must 
be fairly realistic because the scaling factors that give 
the best fit to the spectra at the small and large accep- 
tance angles are essentially the same. The analysis also 
shows that the adopted structure model must be fairly 
realistic because the best fit to the data is achieved with 
a scaling factor that is very close to one. 

Having achieved a comprehensive description of the 
measured decay-energy spectra, the structure model was 
applied to calculate the cross sections for the radiative 
neutron capture on 14 C, not only to the ground state but 
also to an excited state of 15 C. The total capture cross 
section predicted by the model is (accidentally) in perfect 
agreement with the recently measured, low-energy cap- 
ture cross section, and it also accounts for the neutron 
capture data at higher energies. That implies that the 
measured decay-energy spectra and the measured neu- 



8 



tron capture cross sections are consistent within a few 
percent, when the decay-energy spectra are analyzed in 
semiclassical, dynamical description. 

Good agreement between the dissociation and neutron 
capture measurements was also recently achieved within 
the CDCC description [2]. The analysis of the dissocia- 
tion data was based on structure models that are slightly 
different from the model used here, in particular with re- 
spect to d-waves. A comparison of the semiclassical and 
CDCC calculations will be made in the near future [23| . 
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VII. APPENDIX: MAC CROSS SECTION 

The Maxwellian Averaged Capture (MAC) cross sec- 
tion at the temperature kT is defined by (see Eq. (5) of 
Ref. (H) 



2 fdE E E 

CTMAC = ~~7= I T7^ Vrc(E) — eX P (- — ), 



J kT 



kT 



kT' 



(Al) 



in terms of the radiative capture cross section a rc (E), 
which can be derived from the dipole strength distribu- 
tion. For example, the cross section for the El radioactive 
neutron capture to the l/2 + ground state of 15 C is 



4(2tt) 3 h 2 (S n + E) 3 dB(El) 



E 



dE 



(A2) 



where dB(El)/dE is the dipole excitation strength, S n 
is the neutron separation energy, /i„ the neutron- 14 C re- 
duced mass, and E is the relative energy of the neutron 



and the 14 C fragment. A more general expression can be 
found, for example, in section 3.2 of Ref. [7j. 

Inserting into Eq. (Al) the cross section for the cap- 
ture to the l/2 + ground state of 15 C calculated in the 
model defined by the parameters in Table I, one obtains 
the following relation between the MAC cross section at 
the temperature kT — 23.3 keV and the capture cross 
section at the energy kT =23.3 keV, 



cmac ~ 1-46 a rc (kT). 



(A3) 



(If the small capture cross section to the 5/2 + bound 
state is also included, one obtain essentially the same 
relation, namely, with the factor 1.460 replaced by 1.462.) 
Daniel Baye derived the following expression [22j 



a rc (E) m C V~E(1 + siE + s 2 E 2 ), 



(AA) 



which gives a very good parametrization of the low- 
energy cross section for the radiative neutron capture 
from a p-wave to an s-wave bound state. Inserting this 
expression into Eq. (Al) one obtains the expression 



3 

CMAC = 2 a rc 



(kT) 



1 + fsifcT- 



l s 2 (kTf 



1 + SikT + s 2 (kT) 2 



(A5) 



Eight different structure models of 15 C were considered in 
Ref. [22j and they all give the same relation, Eq. (A3), 
between the capture cross section at 23.3 keV and the 
MAC cross section at the temperature kT = 23.3 keV, 
when the parameters s\ and s 2 of the models are inserted 
into Eq. (A5). 

The cross section for the radiative capture to the 1/2+ 
ground state of 15 C predicted by structure model defined 
in Table I can be accurately parametrized by Eq. (A4) 
(at least up to 1 MeV) with C a V~kT = 4.77 fib, si = 
-0.783 MeV" 1 and s 2 = 0.298 MeV" 2 . 
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